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, ^ ' Statistics of the single-particle levels in a deformed Woods-Saxon potential is an- 



alyzed in terms of the Poisson and Wigner nearest-neighbor distributions for several 
deformations and multipolarities of its surface distortions. We found the significant 
differences of all the distributions with a fixed value of the angular momentum pro- 
jection of the particle, more closely to the Wigner distribution, in contrast to the 
full spectra with Poisson-like behavior. Important shell effects are observed in the 
nearest neighbor spacing distributions, the larger the smaller deformations of the 
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^ . The microscopic many-body interaction of particles of the Fermi systems such as heavy nuclei is 

rather complicated. Therefore, several theoretical approaches to the description of the Hamiltonian 
r> ■ which are based on the statistical properties of its discrete levels are applied for solutions of the 



realistic problems. For a quantitative measure of the degree of chaoticity of the many-body forces, 
the statistical distributions of the spacing between the nearest neighboring levels were introduced, 
first of all in relation to the so called Random Matrix Theory [lH5|. Integrability (order) of the 
system was associated usually to the Poisson-like exponentially decreasing dependence on the 
spacing variable with a maximum at zero while chaoticity was connected more to the Wigner-like 
behavior with the zero spacing probability at zero but with a maximum at some finite value of 
this variable. 

On the other hand, many dynamical problems, in particular, in nuclear physics can be re- 
duced to the collective motion of independent particles in a mean field with a relatively sharp 
time- dependent edge called usually as the effective surface within the microscopic-macroscopic 
** niagner@kinr.kiev.ua 



app roximation [g, [7|. We may begin with the basic ideas of Swiatecki and his collaborators [6|, |8| 
131] . In recent years it became apparent that the collective nuclear dynamics is very much related 
to the nature of the nucleonic motion. This behavior of the nucleonic dynamics is important in 
physical processes like fission or heavy ion collisions where a great amount of the collective energy 
is dissipated into a chaotic nucleonic motion. We have to mention here also very intensive studies 
of the one-body dissipative phenomena described largely through the macroscopic wall formula 
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(w.f.) for the excitation energy |8l-ll3| and also quantum results J9|, [ll|, 

w.f. was suggested originally in Ref. [Sj on the basis of the Thomas- Fermi approach. It was 



14l |. The analytical 



re-derived in many works based on semiclassical and quantum arguments, see Refs ISHlTj for 
instance. However, some problems in the analytical study of a multipolarity dependence of the 
smooth one-body friction and its oscillating corrections as functions of the particle number should 
be still clarified. In particular, we would like to emphasize the importance of the transparent 
classical picture through the Poincare sections and Lyapunov exponents showing the order-chaos 
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transitions [10|, Il8l-l20| and also quantum results for the excitation energy [UJ, |13|, |1J]- Then, the 
peculiarities of the excitation energies for many periods of the oscillations of the classical dynam- 
ics were discussed for several Legendre polynomials, see Ref. [ij]. as the classical measures of 
chaoticity. The shell correction method 7], |2l| was successfully used to describe the shell effects 
in the nuclear deformation energies as functions of the particle numbers. This is important also 
for understanding analytically the origin of the isomers in fission within the periodic orbit theory 
(POT) 22|. We should expect also that the deviations of the level density near the Fermi surface, 
like shell effects, from an averaged constant should influence essentially the nearest neighbor spac- 
ing distribution (NNSD). For a further study of the order-chaos properties of the Fermi systems, 
it might be worth to apply the statistical methods of the description of the single-particle (s.p.) 
levels of a mean-field Hamiltonian within microscopic-macroscopic approaches (see for instance 
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Refs 

The statistics of the spacing between the neighboring levels and their relation to the shell 
effects depending on the specific properties of the s.p. spectra, as well as the multipolarity and 
deformation of the shape surfaces should be expected. The quantitative measure of the order 
(or symmetry) can be the number of the single-valued integrals of motion, except for the energy 
(degree of the degeneracy of the system, see also Refs [2^, |27H29|). If the energy is the only 
one single- valued integral of motion one has the completely chaotic system 30| ■ For the case of 
any such additional integral of motion, say the angular momentum projection for the azimuthal 
symmetry, one finds the symmetry enhancement that is important for calculations of the level 
density as the basic s.p. characteristics. 



Our purpose now is to look at the order-chaos properties of the s.p. levels in terms of the 
Poisson and Wigner distributions with focus to their dependence on the multipolarities, equilib- 
rium deformations and shell effects in relation to the integrability of the Hamiltonian through 
the comparison between the spectra with the fixed angular momenta of particles and full for the 
Woods-Saxon potential. 

II. SPECTRA AND LEVEL DENSITIES 

We are going now to study the statistical properties of the s.p. spectra of the eigenvalue 
problem, 

H<j), = sA. H = T + V, (1) 

where H is a. static mean-field Hamiltonian with the operator of the kinetic energy T and deformed 
axially-symmetric Woods-Saxon (WS) potential, 

V^Vwsir) = -— J"' (2) 

1 + exp |[r — R[0)\ /a} 



mu9 



r, 6, if are the spherical coordinates of the vector r. Following Refs [9|, llO|, ll2|, ll3||, the shape of the 
WS-potential surface is defined by the effective radius R{6) given by: 



m = ^ 






(0) 



(3) 



Here, A is a normalization factor ensuring volume conservation, and a^ stands for keeping a 
position of the center of mass for odd multipolarities and Rq is the radius of the equivalent 



sphere, Y^q(9) = \/{2n + 1)/Atc P„ (cos 6) are the spherical functions. P„ (cos 6*) are the Legendre 
polynomials and a is the deformation parameter independent on time. For diagonalization of the 
Hamiltonian with the WS potential ([2]), the expansion over a basis of the deformed harmonic 



oscillator is used as shown in Ref. 13 1. 



Figs [U and |2] show two examples for the full spectra of the s.p. energies Ei and for the fixed 
angular momentum projection m = versus the deformation parameter a for the P2 and P5 
shapes, respectively. The spectra for P3 and P4 are very similar to the P5 case and therefore, 
they are not shown. As seen from Fig. [T|, there are clear shell effects in the full spectra at small 
deformations, approximately at a; < 0.1 for all multipolarities from the P2 shape to the P5 one. 
With increasing deformation a, the shell gaps become less pronounced and slowly changed in 
the region a ~ 0.1 — 0.4 for all these multipolarities. Much more differences can be found in 
comparison of Fig. [1] for full spectra and Fig. [2] for m = levels only. The shell effects are seen 



here too but much less pronounced. The spectrum of levels with m = becomes more uniform 
with increasing multipolarity n. 

The key quantity for calculations of the NNSD P{S) is the level density g{6), see Appendix 
A and Refs [l|-l5[. For these calculations one may apply the Strutinsky shell correction method 
writing 

g^{e) = gie) + 5g^ie). (4) 

The smooth part g{e) is defined by the Strutinsky smoothing procedure [7|, |2l|. The so called 
plateau condition (stability of values of the smooth level density g as function of the averaging 
parameters: Gaussian width F and the degree of the correction polynomial M takes place at 
F = 20 — 40 MeV and M = 4 — 8). Figs [3] and H] for full spectra and for fixed angular momentum 
projection m = show the typical examples of the level densities (smooth component and the 
total density with the oscillating part) for the same degree of the Legendre polynomials n = 2 
and 5 at deformations a = 0.1 and 0.4, in correspondence with spectra presented in Figs [1] and [2], 
respectively. In Fig. H] for the case of the specific m = levels, one has somewhat larger Gaussian 
width parameters of the smooth level density g in the total density (Eq. (jl])) than those for the 
full spectra in Fig. [3l As seen clearly from Figs |3] and H] the smooth level density for the 171 = 
spectra is more flat (besides of relatively small remaining oscillations because of much less levels 
with the fixed m = 0), as compared to the full-spectra results. This more fiat behavior for the 
fixed m value is due to the loss of the symmetry. Therefore one expects the system with the fixed 
angular momentum m to be more chaotic. The differences between Fig. [3] and Fig. H] are quite 
remarkable. On the other hand differences between pictures within Fig. [3] or Fig. H] are less notable 
and as one can see shell effects are still remaining at bigger deformations and multipolarities. 

III. NEAREST-NEIGHBOR SPACING DISTRIBUTIONS 



Following the review paper [3(], the distribution V{S) for the probability o_f finding the spacing 



5* between the nearest neighboring levels is given by (see also Refs 
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and Appendix A) 



V{S) = g{S) exp ( - / g{x) dx 



/K. (5) 



The key quantity g{S) can be considered as the density of the s.p. levels counted from a given 
energy, say, the Fermi energy Ep. D is a mean uniform distance between neighboring levels so 
that 1/D is the mean density of levels. K is the normalization factor for large enough maximal 



value of S", Smax, 



^ = / dxg{x) expi- / giy)dy\ /D. (6) 



This normalization factor K can be found from the normahzation conditions: 

/d.Pw^/d..n.)^i. (7) 

[Notice that for convenience we introduced the dimensionless probabihty V in contrast to that of 
Ref. J3| denoted as P{S), see Eq. (1.3) there.] 

The Poisson law follows if we take constant for the level density, g{S) = 1/D, in Eq. (jS]), 

V{S) = exp i-S/D) . (8) 

Wigner's law follows from the assumption of the linear level density, proportional to S, 

V{S) = {7x8/20) exp {-tiS^/AD^) . (9) 

Both distributions are normalized to one for large enough maximal value of 5, Smax = cxd to satisfy 
Eq. (CD. 

The level density in fact is not a constant or oc S. The combination of the Poisson and Wigner 



distributions was suggested in Ref. [3l| by introducing one parameter. For our purpose to keep a 
link with the properties of the level density, like smooth and shell components [7|, it is convenient 
to define the probability V{S) (Eq. ([5])) for a general linear level-density function through two 
parameters A and i3, 

g{S) = {A + BS/D)/D. (10) 

Substituting Eq. ( ITOl) into the general formula (Eq. ([5])) one obtains explicitly the analytical result 
in terms of the standard error functions, erf(z) = 2 J^ dx exp(— x^)/y^, 

V(S) = (1 + B(,IA) exp (-Be 12 - Al^) / |N„ + B Hi/^] , (11) 



/■■^j, ( B,j ,\ fir (A^\ JA + BC 

Ko = I df exp [--e - A^) = ^^ exp (^ j erf [-^ 



exp( -^C^-AC] +A^o 



(12) 
where ^ = S/ D, C = S„iax/ D is the maximal value of ^. For large C — > oo one has simply 



No -^ \A72^ exp {A^/2B'^) and Ni -^ -A^q/B. Taking the limits ^ -> 1, -B ^ and ^ ^ 0, 
i3 — )■ 1 in ( [TT|) one simply finds exactly the standard Poisson (Eq. (j8])) and Wigner (Eq. (jQ])) 
distributions. In this way the constants A and B are measures of the probability to have Poisson 
and Wigner distributions (Eq. ( ITOj) ). 



IV. NUMERICAL RESULTS 

Figs [5] and [6] show the corresponding NNSD P{S) (Eq. ([5])). Again, in accordance with spectra 
(see Figs [H E]) and level-density calculations in Figs |3] and HI the dramatic changes are observed 
between Fig. |6]for the NNSD P{S) with the 171 = and Fig. |5]for those of the full spectra ones. 
Results presented by heavy dots in Fig. |5] look more close to the Poisson distribution and those 
in Fig. |6] are more close to the Wigner one. 

There are a large difference in numbers A and B which measure the closeness of the distributions 
P{S) for the neighboring levels spacing to the standard ones, Poisson (1,0) and Wigner (0,1). 
However, in Fig. [6] all distributions are more close to the Wigner in shape having a maximum 
between zero and large compared to D value Smax with respect to D [ Smax = CD, see immediately 
after Eq. ( IT2l) ] than monotonous exponential-like decrease similar to the Poisson behavior in Fig. 
[5l Notice that we have more pronounced Wigner-like distribution with increasing multipolarity 
n and deformation a in Fig. |6l especially remarkable at P5 surface distortions and large enough 
deformation a = 0.4, see last plot (d) in Fig. |6l Including all the angular momentum projections 
m for all desired multipolarities and deformations one has clearly Poisson-like behavior though 
they differ essentially in numbers A, B from the standard ones (1,0), see Fig. [6l 



Q 



he reason for this can be understood looking at the Poincare sections shown in Fig. [7] 
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18l |. The upper row is related to a small deformation and lower row corresponds to a large 
deformation. The projection of the angular momentum is ?7i = in all pictures of Fig. [71 Difference 
is remarkable for the integrable spheroidal cavity and other non-integrable (in the plane of the 
symmetry axis) shapes. As seen from comparison of upper and lower plot lines, with increasing 
deformation a and multipolarity n we find more chaotic behavior and we should expect therefore 
the NNSD closer to the Wigner distribution ([9]). This is in agreement with the NNSD calculations 
for the fixed m = 0, see Fig. [6l Notice that similar properties of the NNSD for other potentials 



and constraints were discussed in Refs 
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-|26|. 



The difference between the NNSD calculations with the realistic level densities by the Strutinsky 
shell-correction method (see Eq. (jl])) for the considered WS potential and those with their idealistic 
linear behavior (Eq. (TTO!) ) can be studied in terms of the general formula P{S) (Eq. ([5])). In 
particular, the shell effects related to the inhomogeneity of the s.p. levels near the Fermi surface 
for all desired multipolarities and deformations are found to be significant, also in relation to the 
fixed quantum number m. 

Figs [HI and [9] show the results of these calculations corresponding to Figs [5] and El Notice that 
in the case of the full spectra, see Fig. [HI one has Poisson-like distributions corresponding to the 
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smooth density (dashed) with a similar behavior as NNSD shown by heavy dots, in contrast to 
Fig. [9] where we find rather big differences between these curves. The shell effects are measured by 
the differences between the solid curve related to the total level density with the shell components 
and the dashed one for the smooth level density of Figs [3] (all m) and H] (with m = 0), see 
correspondingly Figs |8] and [91 With increasing deformations one has slightly decreasing the shell 
effects, in contrast to the multipolarity dependence for which there is almost no change of the 
shell effects at the same deformations. 

V. CONCLUSIONS 

We studied the statistics of the neighboring s.p. levels in the WS potential for several typical 
multipolarities and deformations of the surface shapes and deformations, as compared with the 
standard Poisson and Wigner distributions P{S). For the sake of comparison, we derived ana- 
lytically the combine asymptotic Poisson- Wigner distribution P{S) related to the general linear 
dependence of the corresponding s.p. level density. We found the significant differences between 
distributions for a fixed value of the angular momentum projection m of the particle and those 
accounting all possible values of m. For the case of the fixed m = we obtained distributions P{S) 
more close to the Wigner shape with the maximum between 5* = and a maximal large value of S, 
the more pronounced the larger multipolarity and deformation of the potential surface. We found 
also that the full spectra distributions P{S) look Poisson-like in a sense that they have maximum 
at 5* = and almost exponential decrease as a function of the energy near the Fermi surface. Our 
results clarify the widely extended opinion of the relation of the distributions (Poisson or Wigner) 
to the integrability of the problem (the integrable or chaotic one). All considered potentials are 
axially-symmetric but they are the same non-integrable ones in the plane of the symmetry axis. 
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28l /C, i.e. the number of the 



However, the degree of the symmetry (classical degeneracy [22|, 
single-valued integrals of motion besides of the energy) for the case of the full spectra. Figs [H [HI 
[5l ( /C = 1, a mixed system) is higher than for the fixed angular momentum m {K, = like for 
the completely chaotic system). Notice that integrability is not only one criterium of chaoticity. 
The measure of the differences of the distributions P{S) between Poisson and Wigner standard 
ones depends also on the properties of the energy dependence of the level density (from constant 
to proportional-to-energy dependence). From comparison between the general distribution P{S) 
related to the smooth level density obtained by the Strutinsky shell correction method and the 
statistics of the neighboring s.p. levels one finds that all of them are more close to the Poisson-like 
behavior. This shows that the energy dependence of the smooth level density differs much from 



the linear functions. We obtained also large shell effects in the distributions P{S) in nice agree- 
ment with those in the key quantity in this analysis, - level density dependencies on the energy 
near the Fermi surface. 

As to perspectives, it might be necessary to use the combined microscopic-macroscopic ap- 
proaches p, l7| to clear up the results more systematically and analytically. Our quantum results 
can be interesting for understanding the one-body dissipation at slow and faster collective dynam- 
ics with different shapes like the ones met in nuclear fission and heavy-ion collisions. 
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Appendix A: A derivation of the NNSD 

We introduce first the level density, g{E), as the number of the levels dA^ in the energy interval 
[E,E + dE] divided by the energy interval, g{E) = dN/dE. With the help of this quantity 
one can derive the NNSD P{S) as the probability density versus the spacing S between the 
nearest neighboring levels. Specifying P{S) to the problem with the known s.p. spectra of 
the Hamiltonian, one can split the energy interval AE under the investigation into many small 
(equivalent for simplicity) parts AS ^ AE. Each of AS nevertheless contain many energy levels, 
AS* ^ D. Then, we find the number of the levels which occur inside of the small interval AS*. 
Normalizing these numbers by the total number of the levels inside the total energy interval AE 
one obtains the distribution which we shall call as the probability density P{S). Notice that the 
result of this calculation depends on the energy length of the selected AS*. In our calculations, we 
select AS" by the condition of a sufficient smoothness of the distribution P{S). Such procedure 
is often used for the statistical treatment of the experimentally obtained spectrum with the fixed 
quantum numbers like the angular momentum, parity and so on |3|. 

Following mainly Ref. [5[, let us calculate first the intermediate quantity f{S) as the probability 
that there is no energy level in the energy interval [E, E + S]. According to the general definition 
of the level density mentioned above, g{S)dS can be considered as the probability that there is 
one energy level in [E + S,E + S + dS]. Then, 

fiS + dS) = fiS)il-giS)dS) , (Al) 

which leads to the differential equation for f{S), 

d/ = -g{S)dSf{S). (A2) 



Solving this equation one gets 



fiS)=Cexp[- I g{x)dx] . (A3) 



Let P(5')d5' denote the probabihty that the next energy level is in [E + S,E + S + dS], 

P{S)dS = f{S)g{S)dS . (A4) 

Then, substituting Eq. (1A3I) into Eq. (1A4I) one finally arrives at the general distribution: 

PiS) = CgiS) exp f- J giS')dS') . (A5) 

The boundary conditions in solving the differential equation (IA2I) accounts for the meaning of the 
NNSD P{S) and its argument as the spacing between the nearest neighbor levels as shown in the 
integration limit in Eq. (IA5I) . The constant C is determined from the normalization conditions 
(Eq. dSD). 
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FIG. 1. The s.p. energy levels Si in the WS potential {Vq = 200 MeV, Rq = 6.622 fm, a = 0.1 fm) as 
function of the deformation a for the P2 (left) and P5 (right) shapes (Eq. ([3])). 
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FIG. 2. The same s.p. spectrum of levels as in Fig. [T] but with the projection of the angular momentum 
m = 0. 
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FIG. 3. The level densities gr{S) (Eq. (j3|)) as a function of the energy counted from the Fermi level 
{S = 0) for a given particle number N for spectra of the s.p. levels of Fig. [1] for the P2 (a, b) and P5 (c, d) 
shapes at the small a = 0.1 (left, a,c) and larger a = 0.4 (right, 6, d) deformations; dashed is the smooth 
density g{S) (Eq. Q); solid is the total density gr{S) (Eq. ©) (F = 3 MeV, M = for 6griS)). 
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energies S for the same spectra as in Figs [1] and [3l Solid curve is a standard Poisson distribution (Eq. 
dSl)) and a dotted one is a standard Wigner distribution ([9]). Numbers in brackets {A,B) show A and B 
of Eq. (|11|) . Dashed curve corresponds to a linear approximation to the level density (Eq. (|1U|) 1: other 
notations are the same as in Figs [3] and [H 
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FIG. 6. The same as in Fig. [5] but for the spectra of Figs [21 [H 
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FIG. 7. Poincare sections Vp vs p for spheroid, P2 and P5 shapes at the small deformation a = 0.1 
(upper row) and large deformation a = 0.4 (lower row) for the projections of the angular momentum on 
the symmetry axis m = 0. 
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FIG. 8. The general distribution P{S) (Eq. ([5])) vs the energies S for the same spectra as in Figs [U El 
[5j dashed are the distributions P{S) related to the Strutinsky smooth density g{S) and solid is the total 
level density griS) (Eq. ([4])); dots are the same as in Figs [5] and [H 
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FIG. 9. The same as in Fig. [8] but for spectra of Figs [21 H] and [H 



